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Abstract 

We consider two-dimensional t — J — V model. The wave function of the 
ground state is constructed. We demonstrate that the doping by holes results 



£> ■ in condensation of the spin- waves, destruction of the long-range antiferromag- 

netic order and formation of the gap in the spin-wave spectrum. 
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It is well known that the long-range antiferromagnetic (AF) order in Cooper Oxide super- 
conductors is destroyed under the doping by the holes. In the frameworks of two-dimensional 
t — J model the origin of this instability was realized, (see e.g. Refs. [ffl-fl). The instability 



is due to the strong interaction of spin-waves with mobile holes. However a structure of 
the ground state as well as a spectrum of excitations was not understood. In the present 



work we discuss the ways to stabilize the hole-hole pairing, construct the ground state wave 
function as a condensate of spin-waves and discuss the spectrum of bosonic excitations. 

At half-filling (one hole per site) the t-J model is equivalent to the Heisenberg AF model 
which has long-range AF order in the ground state. We consider the doped system starting 
from this ground state which we denote by |0). In spite of the destruction of the long-range 
AF order it is convenient to use |0) and corresponding quasiparticle excitations as a basis 
set in the problem with doping. 

The effective Hamiltonian for t — J model was derived in works J?] [J 

H eff = e k^L^ka + W q( a l tt q + A^q) + H h,sw + H hh . (1) 
ker q 

It is expressed in terms of usual spin- waves on AF background a q ,/3 q (see e.g. Ref. |T0|), 
and composite hole operators (a = ±1/2). The summations over k and q are restricted 



inside the Brillouin zone of one sublattice where 7 q = ^(cosq x + cosq y ) > 0. The spin- wave 



dispersion is u; q = 2<Jl — 7^ pa \/2|q|, at g «C 1. Let us recall that the parameters of t — J 
model are t and J. We set hereafter J = 1, so all energies are measured in units of J. Single 
hole dispersion has minima at the face of magnetic Brillouin zone ko = (±7r/2, ±7r/2). Near 
these points the dispersion can be presented in a usual quadratic form e p ~ hPiPi + \P2P2 
(P2 *C Pi), where p = k — k , and P2 is projection along the face of Brillouin zone, pi is 



orthogonal projection of the momentum. For 5 > t > 1 P\ Prf 0.65t (see e.g. Refs. [[ll"HT3" 



Following Refs. [11,13] we will set a = P1/P2 ~ 5 — 7. For the small concentrations 5^1 
under consideration, holes are localized in momentum space in the vicinity of the minima of 
the band and the Fermi surface consists of ellipses. The Fermi energy and Fermi momentum 
of non-interacting holes are: ep = ^ir(PxP2) 1 ^ 2 8, Pf ~ (ixd) 1 / 2 . The Fermi momentum pp is 
measured from the center of the corresponding ellipse. Let us stress that the numerical value 



of ep is very small. For realistic superconductors t/J ~ 3 (see, e.g., Refs. Hl4|jr5| ). Therefore 
at 8 = 0.1 and J = 0.15 eV one gets e F ^ 15 meV ^ 175 K. 

The effective interaction of a composite hole with a spin-wave is of the form (see, e.g. 
Refs. [0,0]) 

9(K q) + + H.C. ) . (2) 

k,q V 7 

For k pa ko and q 1 the vertex is g(k, q) ~ 2 3//4 /gi/ Jq. The component qi is perpendicular 
to the face of Brillouin zone. The coupling constant / was calculated in the Ref. H. For 
t = 3 it is equal / ~ 1.8, for large t the coupling constant / approaches to 2. 

The interaction between the two holes can be caused by the exchange of single spin-wave. 
Alongside with that there is a contact hole- hole interaction which is denoted in ([J) by Hhh- 
It is of the form 

r c i 

^7k!-k 3 + Y(j k i+ k 3 +7k 2 +k 4 ) ^ST^Il ^2^1^12,34- (3) 

An expressions for the coefficients A and C as a functions of t are presented in the works 



ihh ~ o 

1,2,3,4 



It was demonstrated in the works |16],[T7| that the spin-wave exchange results in very 
strong pairing between the holes. The pairing is strongest in d- and g-waves where the 
superconducting gap A is of the order of Fermi energy ep. The approach works [0,0 is 
based on the observation that at the typical momentum q ~ pp the spin- wave frequency is 
much larger then Fermi energy: uj q ~ pp ~ (ix5) l l 2 ep = \{PiP2) l ^ 2 5. This is why we 



can calculate pairing using unrenormalized spin-waves. Now we are going to discuss what 
happens with spin- waves at q < pp. Let us consider renormalized spin- wave Green function. 

G(uj,q) = — 1 — -. (4) 

uj 2 - uj\ + P(u, q) 

For stability of the system the condition 

wj > p(P, q) (5) 

should be fulfilled. Otherwise the Green function (|j) would possesses a pole with imaginary 
uj. The diagrams for polarization operator are presented in Fig.l. Chaining is due to the 
contact interaction (|3|). In t — J model the constant C in contact interaction is small at any 
t, and the constant A vanishes exactly near the point which we are interested in: t « 3. 
Therefore the contact interaction is small and polarization operator is given by first diagram 
only. If we assume the holes be a normal Fermi liquid the calculation P(uj, q) is very simple. 
For q <C p F one gets 

mq)^i? 2 (6) 

7TVP1P2 

After the substitution of the values of / and /?i , /? 2 presented above, we see that condition 
([5D is violated, and P is about 2.5 times larger than uj 2 . This indicates the instability of 
normal Fermi liquid ground state What happens if we take into account the hole- hole 

pairing? In this case there is no simple analytical expression for the polarization operator. 
Our numerical computations show that pairing practically does not influence P(0, q). With 
pairing calculated in the works ]TJ],|l7|] -P(0, q) is only by 8% smaller then the value given by 
(^). Let us note that it is quite strong pairing: the maximal value of superconducting gap on 
the Fermi surface is A ~ 0.7ep, and this gives reasonable values of critical temperature. Even 
if we enhance pairing by hands up to the value A ~ 1.3ep, P(0, q) is only by 14% smaller 
then the value given by (||). Thus even with pairing we have instability of ground state. Let 
us stress that this conclusion is practically independent of the mechanism of pairing. 

One can believe that nonlinear interaction of spin-waves which is not included into ef- 
fective Hamiltonian ([I]) could stabilize the system. However the maximal magnitude of 
imaginary uj is at q ~ pp. Here the spin- wave residue of unstable mode in Green function 
(fD is very small. It means that actually the system is unstable with respect to spin-sound 
in hole Fermi liquid. The picture is as follows: The Green function (SI) has two collective 
poles: l)the upper pole which originates from initial spin-wave, 2)the lower pole which orig- 
inates from spin-sound. Due to the interaction the spin sound is repulsed down from initial 



spin-wave. The repulsion is so strong that it acquires imaginary frequency. The nonlinear 
interaction of spin-waves practically does not influence the spin sound and therefore can not 
eliminate the instability. We do not see any possibility to stabilize t — J model without 
contact interaction. The true ground state in this case is probably some spiral phase. 

Let us introduce now the hole-hole repulsion V at nearest sites and consider t — J — V 
model. In the effective Hamiltonian ([!]) only contact interaction ([3]) is changed: we should 
consider the constant A as independent parameter of the model (A « V). Now the chain in 
Fig.l becomes essential, and simple estimations show that value A <~ 1.0 — 1.3 is enough to 
eliminate the instability. If J = 0.15eV it means that A ~ 0.2eV. It is hardly believed that 
in realistic systems there is no such a small Coulomb repulsion between the holes at nearest 
sites. The repulsion is probably even larger, but in this case, for the chain Fig.l one has to 
use hole-hole scattering amplitude instead of simple matrix element Phh- Calculation of this 
amplitude is in the progress. 

Consider now the problem of long-range AF order in t — J — V model with V big enough, 
so that the paired hole Fermi liquid is stable. For this question it is convenient to use 
Hamiltonian technique instead of Feynman one because we need explicit construction of 
ground state wave function. The polarization operator P(u, q) in (|j) corresponds to nor- 
malization 2u q spin-waves in the volume. For Hamiltonian approach let us introduce the 
polarization operator H(u, q) corresponding to normalization one spin-wave in the volume: 
P(uj, q) = 2u q Il(uJ, q). The wave function of renormalized spin- wave corresponding to Green 
function (|]) is a combination of a q and /3_ q . To find this wave function write down the 
effective spin-wave Hamiltonian. 

H sw = E(K- n ( w »q))K»q + ^q) +n(a;,q)(a q /3_ q + at/3t q )\ (7) 

q V ' 

The term proportional to c<j q comes from "bare" Hamiltonian ([l]). First "n term" comes from 
diagram Fig. 2a where one spin- wave is annihilated and the other is created. (For simplicity 
we do not present the chain with hole-hole contact rescattering) Second "n term" comes 
from diagrams Fig.2bc where two spin-waves are annihilated or created. Let us note that 
spin- waves have definite values of S z : a q has S z — — 1 and /3_ q has S z = +1. Therefore they 
can appear only in combinations presented in (|7|). One can certainly prove this explicitly 
using the vertex @) and calculating the polarization operator. In the second "n term" the 
spin-waves have the opposite momenta. The vertex (0) is proportional to the momentum. 
Just due to this reason the second "n term" has different sign in comparison with first one. 



Dioganalization of Hamiltonian (0) by Bogoliubov transformation gives the spectrum of Bose 
excitations in the system 

nl = K - n) 2 - n 2 = u\ - 2^ q n(^ q , q). (8) 

This is exactly the equation for the poles of Green function @). To find new ground state we 
have dioganalize (0) at u — 0. As usually for Bogoliubov transformation this ground state 
is of the form 

\gs) oc expf^c q a q /3l q J|0). (9) 

^ q ' 

This is exactly the condensate of spin-waves. 

We started from Neel ground state |0) with two sublattices v4-up and 5-down. The 
difference in magnetization of two sublattices is of the form (see e.g. Ref. [|18|1 ) 

l(S%-S z B ) = l-f -2j2 — (« q a q + (3% - 7q(«q/3- q + a\f3lS) , (10) 

l q ^ q V J 

where 1 — fo ~ 0.303. Using parameters of transformation dioganalizing ([?]) one can easily 
calculate renormalized magnetization 

W^- 5 «W = l-/o-2/(^-^)^, (ID 

where f2 0q = y^ q — P(0, q) At small q due to the superconducting gap f2o q ~ fi q . At 
q ^> Pf the polarization operator vanishes and fioq ~~ * ^q- Therefore the integral in Eq.(|TT|) 
converges at q ~ pf, and we get the estimation for 5S Z coursed by spin- wave condensation 




Here v = is unrenormalized spin-wave velocity and v is renormalize that which follows 
from Eq.(||). If \SS Z \ = 0.303 the magnetization vanishes and one should conclude that 
the long range AF order is destroyed. Due to estimation ( |12D for v/v « 4 it happens at 
5 = 5 C ~ 0.04. Note that the considered effect of enhancement of spin quantum fluctuations 
due to the polarization of fermionic subsystem is similar to the well known Casimir effect in 
Quantum Electrodynamics. 

What happens if 5 > 5 C and magnetization calculated using formula (p4]) becomes neg- 
ative? It means that there are a lot of spin-waves in condensate and we have to take into 
account their nonlinear interaction. We can not do it exactly. Fortunately there is a simple 



approximate way suggested by Takahashi in the work on Heisenberg model at nonzero tem- 
perature ||18|| . Following Takahashi let us impose the condition that sublattice magnetization 
vanishes. 

bs\\(S* A -S* B )\g 8 ) = 0. (13) 

To find the ground state with this condition we have to dioganalize 

H„ = H sw -±v 2 (S A -S B ), (14) 

where H sw is given by © and {S% - S Z B ) by (|10D . Simple calculat ion shows that instead of 
(Bl) we get a spectrum of excitations with a gap 



K = M + v 2 - ( 15 ) 



2 



The average value of magnetization is given by the formula ( p]) with = y^o q 
instead of fioq- We have to find the gap v substituting this formula into condition (|i~3|). Let 
us stress that this condition reflects strong nonlinearity of theory. In essence it is effective 



cutoff of unphysical states in Dyson- Maleev approach (see discussion in the work |19 
Certainly the suggested solution is not exact. It is kind of variational approach 
From the Eqs. (|ITD , (|T3|) , (|15|) we conclude 



v oc 



(V5-JQ. (16) 



However for detailed calculations of the spin-wave gap we have to take into account not only 
the "Casimir" contribution (|i~Tl) , (|i~2|) into spin quantum fluctuation, but also the contribution 
which is due to the spin-wave exchange in hole-hole pairing. Such detailed calculation will 
be presented elsewhere. 

In this work we have introduced the short range hole-hole repulsion V. One can prove 
that very small value of V practically destroys the d-wave hole-hole pairing. However it does 
not influence the g-wave pairing which has the same long-range behaviour as d-wave, but 



quite different short-range one |j|[17|- Therefore the presented scenario favours the g-wave 
pairing. 
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FIGURE CAPTIONS 

FIG. 1. Spin- wave polarization operator. 



Fig. 2. Spin-wave polarization operator in Schrodinger representation: 
a) One spin- wave is annihilated and the other is created. 
b,c)Two spin- waves are annihilated or created. 



